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We show that optical tachyonic dispersion corresponding to superluminal (faster than-light) 
group velocities characterizes parametrically amplifying media. The turn-on of parametric ampli- 
fication in finite media, followed by illumination by spectrally narrow probe wavepackets, can give 
rise to transient tachyonic wavepackets. In the stable (sub-threshold) operating regime of an opti- 
cal phase conjugator it is possible to transmit probe pulses with a superluminally advanced peak, 
whereas conjugate reflection is always subluminal. In the unstable (above-threshold) regime, su- 
perluminal response occurs both in reflection and in transmission, at times preceding the onset of 
exponential growth due to the instability. Remarkably, the quantum information transmitted by 
probe or conjugate pulses, albeit causal, is confined to times corresponding to superluminal veloci- 
ties. These phenomena are explicitly analyzed for four-wave mixing, stimulated Raman scattering 
and parametric downconversion. 



I. INTRODUCTION 



A variety of mechanisms is now known to give rise to superluminal (faster-than-c) group velocities, which express 
the peak advancement of electromagnetic pulses reshaped by material media: 

1. Near-resonant absorption anomalous dispersion in the linear regime of an absorbing medium forms the 
basis for this superluminal reshaping mechanism. 

2. Reduced transmission or evanescent wave formation (tunneling) in passive dielectric structures |^-|^: this re- 
shaping mechanism has been attributed to interference between multiply-reflected propagating pulse components 
in the structure ||^. 

3. Soliton propagation in dissipative nonlinear structures superluminal group velocities can occur in such sys- 
tems via nonlinear three- wave exchanges, as in stimulated Brillouin backscattering in the presence of dissipation. 
They also occur in a nonlinear laser amplifier Q . 

4. Pulse propagation in transparent (non-resonant) amplifying media ||9|,|l0|. Superluminal pulse reshaping in this 
regime has been attributed to either the dispersion |^ or the boundary reflections [ pO[ of the amplifying medium. 

5. Tachyonic dispersion in inverted two-level media: the dispersion in such inverted collective systems is analogous 
to the tachyonic dispersion exhibited by a Klein-Gordon particle with imaginary mass pT] |. Consequently, it 
has been suggested |12[ that probe pulses in such media can exhibit superluminal group velocities provided they 
are spectrally narrow. Gain and loss have been assumed to be detrimental for such reshaping. We note that 
Ref. ]l2| describes an infinite medium and boundary effects on the reshaping have not been considered. 

While Ref. [p^ suggests that tachyonic semiclassical propagation is measurable in optics, it provokes several im- 
portant questions: (i) What other optical processes can give rise to tachyonic behavior, and under what conditions? 
(ii) What is the fundamental origin of tachyonic dispersion, and does it require optical coherence^ This question 
is prompted by our earlier finding that coherence is the key to superluminal group velocity of evanescent EM 
wavepackets (photon tunneling) in linear dielectric media (Hi) Most intriguingly, is it possible to generate 

tachyon quanta and thereby transmit information without causality violation? The subtlety of this question is under- 
scored by the following consideration: If a single-photon input could be transformed by a medium of length L into 
a tachyonic quantum, whose detection probability is localized at times shorter than L/c, then such detection would 
provide information superluminally (!). 

Our search for answers to these questions has resulted in the present theory p^ ], which reveals the conditions for 
the observability of EM wavepackets with tachyonic features in any parametrically amplifying (PA) medium. These 
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conditions are explicitly analyzed for processes such as stimulated Raman scattering (SRS), parametric downconversion 
(PDC) and optical phase conjugation (OPC). Two remarkable findings emerge from the present analysis, which 
employs a new measure of wavepacket reshaping, namely, its time-dependent quantum information (change in Shannon 
entropy): (1) Tachyonic wavepackets transmit much less quantum information per photon than the probe wavepackets, 
and their temporal resolution is much worse than L/c, in accordance with causality. However, their information 
content is confined to superluminally short times ^ well below L/c. This finding may lead to the application of 
tachyonic propagation in communications. (2) Unlike evanescent EM wavepackets in linear dielectric media, tachyonic 
wavepackets are insensitive to decoherence. 



II. TACHYONIC DISPERSION IN INFINITE PA MEDIA 



Our analysis will disclose the role of spatial and temporal boundary conditions on tachyonic propagation in PA 
media. To this end, let us first review the properties of an infinite, one dimensional (ID) PA medium: 
(a) The simplest model consists of two parametrically-coupled fields at different (nondegenerate) frequency bands, 
denoted by mode operators {aj:} and {b^} . If their coupling is effected by a classical pump field at frequency ujq, 
then their interaction Hamiltonian is energy nonconserving , 

Hint = -j?i^(gfcfc'afc6fc'e*"''* - h.c). 

kk' 

It describes two-photon creation or annihilation in correlated ("signal" and "idler") modes, at the expense of the 
classical (undepleted) pump field, whose amplitude determines the coupling gkk'- When the following conditions 
hold: («) the fc-dependence (dispersion) of gkk' is negligible, gkk' 9', {H) phase matching prevails, and {Hi) the 
uncoupled frequencies are LOa(b)k — ± c(fc — fco), we can diagonalize this hamiltonian in the interaction picture 
(upon transforming away the ujo- and fco-dependence) by the Bogoliubov transformation (BT) 

Dk = Wfeflfe + Vkbtk, [Dk, H] = huJkDk, (1) 
Ufc=sinhsfc, Wfe = e*'^ coshsfc, tanh^ Sfc = l^l^/jw/t + c/c|^. (2) 

Whereas the significance of Sk as the squeezingparameter is well known [ p^ , much less attention has been paid 
to the dispersion relation obtainable from eqs.(|l|-||) (Fig. 1) 



L.k^±Vicky-\9\' (3) 

which is tachyonic, i.e., isomorphic to that of a Klein-Gordon particle with imaginary mass p^Jl^ , and yields super- 
luminal group velocities \duJk/dk\ > c for spectrally narrow photonic wavepackets. 
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FIG. 1. Dispersion relation (solid line) in a OPC. The dashed (dotted) lines correspond to the dispersion relation for the 
probe (p) and conjugate (c) waves in vacuum and fco = wo/c. 



2 



SQuO VII 



(b) This behavior should be contrasted with that obtained when we replace the parametric-coupling by the energy- 
conserving 

-ffint = ihg'^{akbl - h.c). 

k 

This model pertains to light scattering in passive, non-dissipative media, e.g., periodic dielectric structures with 
two-band dispersion. The diagonalizing BT in this case describes passive light reflection 

Dk = Uhttk + Vkb^k with Uk — sm9k, Vk — e'''^ cos 9 k and ta.v? Ok — \g\'^ /\ijOk + ck\'^ . 
This BT yields the polaritonic dispersion 

Wk = ±^{ckY + (4) 

which is analogous to that of Klein-Gordon particles with real masses ("bradyons") and allows only subluminal group 
velocities \dioldk\ < c. 

A more general model of infinite PA media accounts for both energy conserving and nonconserving exchanges, and 
therefore pertains to a wider range of optical parametric processes 

J2 {f^ckia+^a,k + ^) + hnkib+^b,k + ^) (5) 

j'.J = ±,k>0 

+ j%fe[(a^fe6ife - h.c.) + {ujkbj'k - h.c.)]} 

When this 4-mode coupling Hamiltonian is diagonalized by means of a BT, it yields a 4-branch dispersion, with two 
polaritonic and two tachyonic branches. In the nondispersive limit fi^ ^ fioi this model pertains to SRS with Stokes 
and anti-Stokes scattering corresponding to the energy-conserving and non-conserving terms in eq.(H). The distance 
between the Stokes and anti-Stokes pairs of branches is determined by f^o, whereas each pair is split by g (Fig. 2). 
Strong non-degeneracy of the two pairs of branches fig > g is required to separate between tachyonic and bradyonic 
(polaritonic) types of dispersion (described by eqs.(||) and (||)). Otherwise, hybrid (mixed tachyon-bradyon) branches 
can arise. When strong non-degeneracy exists, the tachyonic branches in Fig. 2 are rotated by 45° relative to those 
corresponding to eq.(||) and are isomorphic to their tachyonic counterparts in the inverted two- level model [^21. 
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k 

FIG. 2. SRS dispersion obtained from a sum of the energy-conserving and nonconserving two-mode hamiltonians discussed 
in text. Upper curves - Stokes dispersion (tachyonic), lower curves - anti-Stokes dispersion (bradyonic or tardyonic). 

The observability of the above tachyonic dispersion poses not only the unrealistic requirement that the PA medium 
be infinite, but also that we detect the field "dressed" by the medium, obeying eqs. (0-||) (or their equivalent for 
eq. (H)), rather than the "bare" field (unperturbed by Hint)- If the PA coupling is turned on at i = throughout a 
very long medium, then a probe wavepacket in a multimode (Gaussian-distributed) state will be reshaped at t > in 
the dressed-field basis to exhibit superluminal peak advancement, without any gross distortion. However, it will be 
accompanied by broadband excitation of quanta arising from vacuum fluctuations within the phase-matched spectral 
band, with relative weights [D^Dk) = cosh^Sfc. In the bare- field basis, the reshaping of the probe wavepacket at 
<: > will be more complicated, and eventually result in its exponential growth. 
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III. QUANTUM DESCRIPTION OF FINITE OPC 



In the rest of the paper, the above reshaping and dispersion wiU be compared with those of probe photons injected 
into finite PA media of length L. We shall adopt a one-dimensional description of propagation, consistently with 
the standard paraxial approximation. Since we wish to explore time-resolved information transfer, the first generic 
regime to be addressed is that of abruptly switched-on pump fields and short (transient) wavepacket transmission 
times L/c'^t <^ TicI, where Trci is the medium relaxation (or dissipation) time. Superluminal features in this regime 
require the simultaneous pumping of the entire medium, otherwise pump retardation dictates the probe propagation 
velocity. Simultaneous pumping of the medium is possible in OPC, where pump pulses can be perpendicular to the 
medium axis z (Fig. 3- inset). 
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FIG. 3. Conjugate (reflected) wavepacket reshaping in OPC: dependence of flux (ric) (in units c/L) a.t z — on time (for 
\g\L = 1.7, probe width (FWHM for Gaussian amplitude) At — 24L/c, probe detuning from PA resonance S = Q.28c/L, 



tp = 50L/c, {Np) = 10 ); 1 - Gaussian probe transmitted in empty space; 2 



corresponding conjugate (ni'^'^''^' ) : 



observe the 



superluminally advanced peak at t ~ 42L/c, followed by exponential growth (instability) at t > 50L/c; 3 - exponentially-growing 
"noise" (ni"^'); 4 - abruptly switched-on probe (onset at lOL/c, delay tp — 2QL/c), 5 - exponentially-growing conjugate 
corresponding to the latter. 



The quantized analysis of OPC following the abrupt turn-on of the pump throughout the medium has been under- 
taken by us for the first time. It amounts to solving the coupled-wave equations for the slowly varying probe (p) and 
conjugate (c) field operators E^'^^^{z,t), 
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(+) 



dz 



at 



(6) 



Here the coupling constant g is proportional to the pump intensity, which is assumed to be turned on at t = 0. The 



non-standard source terms, i.e., Ep \z,t = 0) and Ec^'{z, t — 0), account for the initial {t = 0) vacuum fluctuations 
of the probe and conjugate fields, respectively, at any point z of the medium [Q < z < L). Additional sources of 
vacuum fluctuations are the input fields Ep{z = Q,t) and ErXz = L,t) at the forward and rear faces of the medium, 
respectively, which are taken to be in the vacuum state. The nonvanishing correlation functions of the source terms 
can be concisely expressed as 



{Z, t)Elj>^ (Z', t')>vac - A5{z - z' T cit - t')), 



(7) 



where the upper (lower) sign refers to the probe (conjugate) field, and A = TiL0Qk'^Q,/2'iT , VL being the solid angle in 
which the output is detected. 

We have solved Eqs. (||), taking account of Eq. (^), to allow for spontaneous parametric generation of both probe 
and conjugate fields from vacuum fluctuations, along with response to an injected probe field £pfp{t). Specifically 
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will be a bell-shaped function with the maximum at satisfying |/p(<)| < |/p(io)| = 1- The average conjugate- 
reflected flux (photons/sec) through the cross-section S and solid angle at z = is the sum of the two above 
components, (ric) = {uc^^) + (ric'"'^'*^^), 



{n[^^Hz = 0,t)) = K^a- 



Gl(C,T)dC+ / Gl{0,T')dT' 



Gp(0,T-r')/p(r')dr' 



(8a) 
(8b) 



where k = \g\L is the dimensionless coupling strength, a — HS/Xq, Xq is the pump wavelength, C, = z/L, t = ct/L, 
and (up) is the maximum probe flux, (rip) = (Np)/ \ fp{t)\'^dt, Np being the number of photons in the probe pulse 
which pass through S into the solid angle Q. Solution (0) is expressed in terms of the Green function (App.A) 



G'p(C,r)=.-^ 




(9) 



Here = 2j + ( and = 2(j + !) — (■ are, respectively, the forward- and backward-propagation retardation times 
after j round trips, hjO stands for the modified Bessel function of 2_;/-th order, and 9{) for the Heaviside step function. 

The first term in {n''c^\t)), Eq.(|8^), is the response to the initial (t = 0) vacuum fluctuations of the probe. The 
second term therein represents the response to the probe vacuum fluctuations (taking them as J-correlated) at the 



input (forward) face z = 0. The phase-conjugate response {n, 



(resp) \ 



(Eq.(pb|)) to the incident probe photon flux is 



identical with the classical response to the same probe intensity and pulse shape fp{t). The Green function in (||) 
ensures causality by the step functions 9(t — Tp). These step functions imply that the j = term of the response 
allows superluminal reshapingonly via the forward tail of the probe wavepacket extending throughout the medium 
at t = 0. The j terms in (|^) account for multiple reflections at the boundaries. 

Above the threshold for parametric instability, the spontaneous parametric radiation {ni^^\t)) (triggered by the 
probe vacuum fluctuations) grows at i > nearly exponentially. Nevertheless, there can be quasi-tachyonic reshaping, 



namely, the transformation of a Gaussian probe wavepacket into a wavepacket 



(reap) 



(t)), which resembles a Gaussian 



with a superluminally advanced peak at times < t < t„. The time t — t^ marks the end of quasi-tachyonic reshaping 
and the start of the unstable classical response, i.e., the onset of exponential growth of (ni^'^''^\t)) , so that at t > 
the output becomes unrelated to the input. This growth is changed at i > TicI, TjcI being the relaxation time (see 
below). This sharp separation of the two evolution stages: < t < the tachyon transient time, and the instability 
growth time < t < t^^i is contingent on the following conditions. We need a narrow spectral width of the incident 
probe, large detuning of the probe central frequency from parametric resonance, and large delay time tp between the 



pump turn-on t 
(Fig. 4). 



and the arrival of the incident-probe peak at 2: = 0, otherwise there is fast onset of the instability 
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FIG. 4. Unstable (above-threshold) regime of the OPC: reflection and transmission of an incoming probe pulse. The OPC 
medium is located between < a; < L (see shaded area). The upper (lower) half of each time plot shows the probe (conjugate) 
pulse \£p{x,t)f (- \£c{x,t)\'^), which is moving to the right (left). Parameters used are So = 0.31 c/L, At = 19.3 L/c, gL = 1.7. 
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IV. QUANTUM INFORMATION TRAVEL IN OPC 



What is the physical significance, if any, of the quasi-tachyonic reshaping at < i < <„? To answer this question, we 
have considered a quantity that has not been used before to study wavepacket reshaping, namely, the time-dependent 
quantum information function AHc(t), associated with the Shannon entropy of the detected flux (ricit)). This function 
reads @ 



AH{t) - (n,) log2 1 + 




(n(^P))log2 1 + 



(10) 



A similar function AHp{t) can be defined for the transmitted probe {np{t)). When the probe wavepacket envelope 
fp(t) has a sharp front, these information functions behave causally, i.e., they grow only at times satisfying retardation. 
However, when fp{t) has an extended forward tail, then these functions exhibit apparent supcrluminal advancement 
compared to the same function for a probe in free space (Fig. 5). The reason for this is that the forward tail of the 



probe is transformed into 



(resp)\ 



right after the pump turn-on at t = 0, hence this information growth is due to 



the nonlocal (simultaneous) pump effect throughout the medium, and no retardation is expected. The information 



increases until 



,(sp)\ 



to {n. 



which grows exponentially and contributes negatively to AH, becomes appreciable, compared 



This causes AHc to be confined to earlier times than the peak of {nc{t))\ 




10 20 30 40 50 



t/(L/c) 

FIG. 5. Information transfer AHc (Eq.([ic|)), same parameters as in Fig. 3, except for (Np) — 100; A - for Gaussian probe 
in empty space AH^ follows probe shape; B - for corresponding conjugate (with "noise"), AH^ is peaked at an earlier time 
t ~ 30I//c, than curve 2 in Fig. 3, but is 10^ times smaller; C - for chopped probe (onset at lOL/c) transmitted through 
empty space; D - for conjugate (with "noise") corresponding to the latter, note that AHc is no longer advanced, but peaked 
at t ~ 20L/c, just as curve G. 



Note, however, that the total amount of information transfer is sharply reduced as compared to that of a probe in 
free space, and does not exceed that of the front tail of the free-space probe, as required by causality. Nevertheless, 
the temporal compression of the information due to the nonlocal character of the response is a potentially useful novel 
feature for communications: it indicates the appropriate time intervals for modulation of the probe envelope, so as to 
optimize the information transfer. As seen from Fig. 5 the Gaussian probe envelope can be switched off well before it 
reaches its peak, since the transferred information AHc is derived only from the forward tail! Likewise, the minimal 
useful interval between successive probe peaks should exceed the temporal width of AHc- 



V. TACHYONIC PROPAGATION IN SRS 



In SRS, by contrast to OPC, the pump cannot be switched on simultaneously throughout the sample rather the 
pump is moving along the sample (in the z direction) with a front velocity c. Therefore the Stokes (signal) field in 



6 



SQuO VII 



SRS or signal in PDC has superluminal features only if the pump front sweeps the sample axis z at a speed exceeding 
c, as allowed by the experimentally feasible "searchlight effect" . 

On the other hand, the SRS regime in which superluminal features are displayed even for a stationary pump is 
that of long times, at which transient effects due to the pump turn-on play no role already, but relaxation effects are 
no longer negligible. The slowly varying Stokes field operator in SRS for a stationary laser pump intensity II, and 
relaxation (dissipation) rate F = l/rrd, is given by the known solution 

-ie-^* I dz'Q\z~ z',0)Io[{az't)^^^]-i f dt' f dz'e-^'' Io[(az't')^/^]F\z' ,t')\ . (11) 
Jo Jo Jo J 

Here a = gT = Ak^Il, F{z, t) is the 5-correlated (in space and time) Langevin operator describing the field dissipative 
fluctuations, and Q{z, 0) is the initial spatially 5-correlated atomic inversion operator. 

The plots (Fig. 6) show the output Stokes (signal) of the average intensity (n^) = {Ei~\z,t)Es^\z,t)) obtainable 
from (^ij) for incident Gaussian input wavepackets with appropriately chosen width, mean frequency and peak delay 
(after the pump turn-on at t — 0), reveal a surprising finding: quasi-tachyonic reshaping obtains for large transit times 
t ^ Tici (if the response contribution to (n^) strongly exceeds the vacuum fluctuations contribution, which comes from 
the last two terms in (^l])). This reshaping is essentially similar to that of the transient {t ^ Trei) regime of OPC 
(Fig. 3). Hence, relaxation (and the resulting destruction of temporal coherence) does not preclude the tachyonic 
reshaping, as opposed to the superluminal reshaping of optical evanescent waves, which is essentially a coherent 
(interference) process Here, by contrast, tachyonic reshaping in the presence of dissipation originates from the t'- 
dependence of the integrand e^^*'/i [{azt')-^/'^]/y/t' in ([Tl|), which is similar to that of a Gaussian (spectrally dispersive) 
function. This conclusion can be generalized to any PA medium: for t ^ Trei, it is possible to observe tachyonic 
dispersion by taking advantage of the high spectral selectivity of the parametric gain coefficient, ~ g/[l-|-Ti-ci(w — wq)], 
which is limited to a band of width F = I/ticI- After the pump is switched on, at times t > Tici, the Fourier- 
transformed response function has the generic form G{S) ~ e i+'^rd-')^ ~ g-V(i-('^r<!i*)^), where S — uj — cjq- The 
nearly-Gaussian shape of the response G{S) allows for tachyonic features, when G{S) is convoluted with a Gaussian 
probe. 




FIG. 6. Transmitted SRS photon flux {(us) /T) as a function of time (Ft): dotted line - exponentially-growing spontaneous 
(vacuum fluctuations) "noise"; dashed line - the probe flux in free space; solid line - induced response with superlumi- 
nally-advanced peak compared to dashed line, (a) Detuning S/F = 10, probe pulse width A/F = 0.5, gL — 25; (b) detuning 
5/r = 15, probe pulse width A/F = 0.25, gL = 12. 



VI. CONCLUSIONS 



Our theory has yielded the following conclusions, which elucidate the nature and observability of tachyon effects 
in quantized PAs: (a) Tachyons are the characteristic quanta generated in any infinite parametrically amplifying 
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(PA) medium, involving processes such as stimulated Raman scattering (SRS), parametric down conversion (PDC) 
and optical phase conjugation (OPC). The basic connection between parametrically amplifying (PA) processes and 
tachyonic dispersion is that their inherent energy-nonconservation entails lack of symmetry under time reversal, along 
with translational invariance (momentum conservation or phase matching). This signifies that the space and time (or 
energy and momentum) axes in PA dispersion arc interchanged as compared to the energy-conserving, timc-rcversible 
polaritonic dispersion, corresponding to the 90° rotation of axes when passing from a bradyonic Lorcntz transformation 
to the tachyonic generalized Lorcntz transformation [ p^ . (b) In a PA medium of length L, only probe wavepackets 
with macroscopic photon numbers and temporal width well above L/ c can be reshaped into observable (transmitted 
or reflected) tachyonic wavepackets (characterized by superluminal peak advancement and little distortion). The 
observability of these tachyonic wavepackets is limited to transient times by the inherent quantum instability of PA 
processes, i.e., the broadband amplification of vacuum fluctuations on the one hand, and, on the other hand, by the 
classical instability (exponential growth) associated with gain, (c) As a result of competition between reshaping and 
instability the information content of tachyonic wavepackets is conflned to times well below L/c, although causality is 
not violated. This novel feature is shown to be advantageous for communications, (d) Tachyonic dispersion of probe 
wavepackets is observable even in the dissipative, incoherent propagation regime, by virtue of the PA gain dispersion 
in the presence of dissipation. 

We acknowledge useful discussions with I. Bialinicki-Birula and Y. Silberberg, and the support of Minerva and 
EU(TMR) grants. 



APPENDIX A: THE GREEN FUNCTION DERIVATION FOR OPC 



In this appendix we briefly outline the two-sided Laplace transform (TSLT) technique introduced by Fisher et al. 
for the derivation of Green function. The TSLT is defined as T(^^)(x,s) = f^^dt F'-*'> {x,t) e'"* . It is only valid 
for functions F(x,t) that diminish faster than exponentially at times t — ^ — oo. Its advantage compared to the the 
usual one-sided Laplace transform is that it also applies to functions that do not vanish at i < 0. The starting point 
of the analysis is to apply the TSLT to the four-wave mixing equations (H) with F(x,t) ~ £p{x,t) or £*{x,t). We 
then obtain the coupled equations 

4z£p{x, s) + ^£p{x, s) + ig£c.*{x, s) = 

(Al) 

■^£c,*ix,s) - ^£c,*ix,s) +ig* £p{x,s) = 0. 



Equations ( Al) are solved together with the Laplace transforms of the boundary conditions fp(0, t) = F{Q, t), where 
F(0, t) is an incident probe pulse at the entry a; = of the PCM medium, and £c{L, t) — 0, so no incoming conjugate 
pulse at the end of the medium. The result is 

£c{x, s) — hr{x, is) .F»(0, s) 

£p{x, s) = ht{x, is) T{0, s), 

with the reflection and transmission amplitudes 

hAx,6)^ , \9WHP{L-x)) 

^sm{f3L) +i(3cos{(3L) ^ ' 

z/3cos(/3(L-x)) + |sin(/3(L-x)) 
^sm(/?L) -I- ipcos{pL) 

where 



f3=-^5^ + \gc\^. (A4) 
c 

The reflected phase-conjugate pulse £c{x, t) and transmitted probe pulse £p{x, t) at position x in the nonlinear 
medium and time t are then obtained by using the inverse Laplace transform. At x = and x = L, respectively, they 
are given by 



1 



7 + ZOO 



^c(0, t) = — ds r,p{is) £p. (0, s) e'' (A5) 

Jj-ioci 



1 



7-4-zoo 



£p{LA)^— dstpp{is)£p{0,s)e'\ (A6) 

Jj-ioo 
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fe(0,i) ^ J dt £;{0,t)H,{t,t) 
£p{L,t) = [ dt £p{0,t)Hp{t,t'), 



where Vcp = hr{x = 0) and tpp = ht{x = L). The choice of contour 7 is in agreement with causahty, which means in 
this case that it is to the right of aU the singularities of rcp{is) (and tpp{ is)) . The sing ularities in the right-half s-plane 
give rise to exponential growth of £c{x,t) and £p{x,t). The integrals (A5) and (A6) can be evaluated by taking the 
contribution of these poles separately and rewriting the remaining integral as a Fourier integral. 

They can also be evaluated in another way, whi ch is especially insightful when regarding non-analytic, chopped 
probe pulses. To that end we rewrite (M) and ( |A6| ) as 

(A7) 
(A8) 

(A9) 
(AlO) 

(All) 



with 



Hpit, t 



1 

1 

27ri 



7+100 



dsr^p{is)e-'^' 



7—100 

7+ioc 



ds tpp(is) e 



-~s{t -t) 



7 — JCJO 



In order to evaluate the integral in, e.g., Hp{t, t ), tpp is first rewritten as the series 



tpp{is) = 2ry^ 



^2n 



-(n+i)r;r 



"0 + v) 



2n+l 



with a = kqc, 77 — \/ — and r ~ 2L/c, the roundtrip time of th e PC M. One can then easily prove that ( All ) is 
uniformly convergent, which allows for term by term integration of (AlO). We use the substitution 



la I u 



An 
U 



with 



't-t- {n+\)T 
t-t' + {n + i)T 



1/2 



(A12) 



(A13) 



in which t + {n + i)r is the retardation time after (n + i) round trips. We then arrive at the following form of the 
conjugate reflected pulse 



£ciL,t) 



dt £;{0,t)Go{t,t) 



00 rt-riT 



-2E 



dt£;{0,t)Gn{t,t), 



(A14) 



with 



and 



Gnit,t) 



—I KqC 

^ 4 



S,r^/2n-2 \g\cy/{t-t'y~n^ 



\g\c^it^t')^-n^ 



n ^2n+2 



(A15) 



\g\c,/{t^t'y-n^T^] ) 



Bn — 



t — t — nr 

t — t' + TIT 



(A16) 
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